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Abstract 

In view of the appreciable semiconducting gap of 0.26 eV observed in recent experiments, epitaxial graphene on 
a SiC substrate seems a promising channel material for FETs. Indeed, it is two-dimensional - and therefore does not 
require prohibitive lithography - and exhibits a wider gap than other alternative options, such as bilayer graphene. 
Here we propose a model and assess the achievable performance of a nanoscale FET based on epitaxial graphene 
on SiC, conducting an exploration of the design parameter space. We show that the current can be modulated by 4 
orders of magnitude; for digital applications an Ion/ loS ratio of 50 and a subthreshold slope of 145 mV/decade can 
be obtained with a supply voltage of 0.25 V. This represents a significant progress towards solid-state integration of 
graphene electronics, but not yet sufficient for digital applications. 

I. Introduction 

Epitaxially grown graphene on SiC provide potential for large scale integration of graphene electronics. The first 
challenge to the use of graphene as a channel material for FETs is to induce a reasonable gap for room temperature 
operation. Recently Zhou et al. [1] have experimentally demonstrated that a graphene layer, epitaxially grown on 
a SiC substrate, can exhibit a gap of about 0.26 eV, measured by angle-resolved photo-emission spectroscopy. The 
gap is probably due to symmetry breaking between the two sublattices forming the graphene crystalline structure, as 
also confirmed by recent density functional calculations [2], [3]. According to the authors of Ref. [1], this method of 
inducing a gap is very easy and reproducible; in addition, the thickness of graphite grown on SiC can be precisely 
controlled to be either single- or multiply layered depending on growth parameters [4]. From a manufacturability 
point of view it is also extremely promising, since it would be highly convenient to prepare an entire substrate of 
graphene on an insulator and then obtain single device and integrated circuit through patterning [5]. 
From its isolation [6], [7], graphene has attracted the attention of the scientific community due to its exceptional 
physical properties, such as an electron mobility exceeding more than 10 times that of silicon wafers [8], and in 
view of its possible applications in transistors [9] and in sensors [10]. To induce a gap in graphene structures, 
several methods have been used: lateral confinement in graphene ribbons [9], [11] carbon nanotubes [12], impurity 
doping [13], or a combination of single and bilayer graphene regions [14], [15], [16]. Unfortunately, they all face 
different problems. 

Carbon nanotubes exhibit large intrinsic contact resistance and are difficult to pattern in a reproducible way; the 
inability to control tube chiraUty, and thus whether or not they are metallic or semiconducting, make solid state 
integration still prohibitive. Graphene nanoribbons [9], [11] allow to obtain a very interesting device behavior [17], 
but require extremely narrow ribbons with single-atom precision, since a difference of only one dimer line in the 
width may yield a quasi-zero gap nanoribbon. Bilayer graphene exhibits a gap in the presence of a perpendicular 
electric field, but the range of applicable bias can only induce a gap of 100 — 150 meV, not sufficient to obtain a 
satisfactory behavior in terms of lon/IoS ratio [15]. Graphene on SiC can is a two-dimensional material, thus does 
not require extremely sophisticated lithography, and provides a higher energy gap: for the sake of comparison, a 
0.26 eV energy gap would require an armchair nanoribbon of width smaller than 3 nm, or nanotubes with diameter 
smaller than 2 nm. 

In this work we present a semi-analytical model of an FET with a channel of epitaxial graphene grown on a SiC 
substrate, where the band structure, the electrostatics, thermionic and band-to-band tunneling cuiTents are carefully 



September 25, 2009 



DRAFT 



2 



accounted for. On the basis of our model, we assess the achievable device performance through an exploration of 
the device parameter space, and gain understanding of the main aspects affecting device operation. 

II. Model 

We adopt the Tight Binding (TB) Hamiltonian for single layer graphene on SiC that was proposed by Zhou et 
al. [1]. The empirical TB valence (— ) and conduction (+) bands of a single epitaxial layer of graphene on SiC, 
read: 

where t is the in-plane hopping term (2.7 eV), m = 0.13 eV is an empirical potential energy shift between the two 
inequivalent graphene sublattices due to interaction with the SiC substrate, and /(k) is the off-diagonal element of 
the considered Hamiltonian [1]. In the six Dirac points of the graphene Brillouin zone, where /(k) is zero, there 
is a finite energy gap Eg — 2m, corresponding to the channel conduction minimum Ecc = m — q(t>ch and the 
channel valence maximum Eye — ~m — q(l)ch, where q is the electron charge and (pch the self-consistent potential 
in the central region of the channel. 

The device under consideration, depicted in \lla), is a transistor with a channel of epitaxial graphene on a SiC 
substrate of thickness tgub = 100 nm, with a top gate separated by a Si02 layer of thickness tox- InEb) we have 
sketched the band edge profiles along the transport direction x, where Eci and Evi respectively represent the 
conduction and valence band edges in the three different regions denoted by i=S, D, C (Source, Drain, Channel). 
Source and drain contacts are n+ doped, with molar fraction an, which translates into an energy difference A 
between the electrochemical potential if^o) and the conduction band edge Ecs (Ecd) at the source (drain) 
contact. The potential is set to zero at the source and to Vds at the drain contact. In the center of the channel 4>ch is 
imposed by vertical electrostatics. We assume, as usual, complete phase randomization along the channel, which is 
particularly important because it allows us to neglect the effect of resonances in the presence of tunneling barriers. 

Exploiting the Gauss theorem we can write the surface charge density in the central part of the channel as 

Q = -Cg (Vg - VpBt - 4'ch) - Csub (Vsub " VpBb " </'c/l ) , (2) 

where Cg = esio^/tox {Csub = ^Sic/tsub) is the capacitance per unit area between the channel and the top gate 
(back gate), Vg (Vgub) is the top gate (back gate) voltage, Vpst (VFBb) is the flat-band voltage of the top gate 
(back gate), which we set to —0.4 eV. 

The transit time of the device in the channel has been estimated as = Qth^c ^ 10"^^ s where Qth and Jth 

''th 

are the thermionic charge and current, respectively, Lc = 20 nm is the channel length. 

In certain spectral regions, for example in the valence band when the device is in the off state, carriers are quasi 
confined by tunneling barriers, and can dwell in the channel for a much longer time and be subject to some degree 
of inelastic relaxation, even if transport in the conduction band is practically ballistic. To consider this effect, we 
have therefore included a degree of inelastic scattering that leads to energy relaxation. 

In steady-state conditions, considering an infinitesimal element of area dkxdky in the wave-vectors space, charge 
distribution in the channel is obtained as a balance between two types of charge exchange processes with the 
contacts: one elastic, and one inelastic. 

We can write the electron charge in the channel as the sum of two contributions: N"^ and Nt- N'^ = f1Lcn2D 
(Nt = flLcn2D) represents the density of forward (backward) going electrons, /|. (kx, ky) denotes the occupation 
factors of forward (+) and backward (— ) states in the channel and {kx, ky) is the 2-dimensional density of 
states in the /c-space. For each contribution we can write a rate equation in steady-state conditions: 

dN"^ _ fe 

-if = 4 - 4 + (1 - Ts) f-vx - (1 - Td) r+vx + '^ Lcn^D = (3) 
at TS 

^ = - J- + J- - (1 - Ts) flvx + (1 - Td) r+Vx + -'- Lcn^D = (4) 

at T£) 

where: ffjs^ -S {E+ — q(f>ch — I^d.s) is the occupation factor at the drain (D) and source (S) contacts, and ^ is 
Fermi-Dirac distribution function. 

Let us focus on eq. [3] (similar considerations can be made for eq. H]): Jg = TsfsVx is the tunneling current 
component injected from source, — TdJ^Vx is instead the drain tunneling current component ejected to the 
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Fig. 1 . a) Schematic picture of a graphene on SiC transistor. The grey line between Si02 and SiC oxide represents the graphene plane acting as 
device channel. Source (S) and drain (D) contacts are also in graphene. b) Protile band structure along the transport direction. The dashed lines 
mark the energy region in which it is possible to have thermionic current Qth), tunneling current from source to channel Oto) and tunneling 
cun'ent from drain to channel Ots)- 



drain, Tg^o are the transmission probabilities from source/drain contacts to channel and is the group velocity. 
(1 — Ts) flvx and (1 — Td) f^v^ are the reflected current components from source and drain barriers, respectively. 
The last term of eqs. [3] and |4] is a thermaUzation process with the source and drain reservoirs, with characteristic 
times Ts and td, respectively. The steady-state and can be obtained by solving eq. [3] and eq. H) 

/+ = ^ 7 — ^7 — , (5) 



,e {^s + ^) (1 - Tn) /I + + 1) {Tu + ^) fh 

f- = ^-^ — , (6) 

where v = i^^" is the inverse of the crossing time Tt- The same reasoning can be applied to derive the hole 
occupation factors in the channel f^.. 

The charge, to be self-consistently solved with eq|2] in order to obtain the channel potential (j)ch, is computed 
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through the integration on the BZ 

where the total current density is expressed as [18] 

*°* - 47r2 



BZ 



BZ 



in 
in 



/_j dkxdky^ 



(7) 



(/+ - /-) rffc^rffcy 



V. 



(8) 



The transmission probability Ts {To) of the interband barrier at source (drain) is zero in the source (drain) 
band gap and 1 when there is no barrier between source (drain) and channel. When a barrier is present Ts is 
computed analytically with the WKB approximation, assuming ky conservation due translational invariance along 
the y direction: 

f /"^'^ _ 1 

(9) 



Ts(£;,A:j,) =exp|-2 j |/m(fc^''='' (x))|da;| , 

where xi and X2 are the classical turning points, and E is the particle kinetic energy. The same approach is repeated 
for Td. 

The potential profile between each contact and the central region of the channel is described by an exponential, 
with characteristic variation length A, obtained from evanescent mode analysis [19]. Assuming tsub ^ tox > ich 
we obtain: 

tch'' 



A 



(10) 



where tch is the effective separation between the interfaces of the Si02 and SiC layers, for which we assume 
tch = 1 nm [9]. 

III. Electrostatics 

From analysis of the electrostatics we can gain a better insight of the device performance limitations. In fact gate 
voltage control upon the channel potential (of which the subthreshold slope 5 is a measure) is strictly limited by 
the quantum capacitance Cq of the channel. Device electrostatics can be schematized as in |2a). The differential 
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Fig. 2. a) Equivalent circuit of device electrostatics, b) Quantum capacitance-voltage characteristics for Vda = 0,0.1,0.25 V, tox = 1 nm, 
taub = too nm and ao = 9.3 X 10"'^ in case of fully ballistic transport. 



capacitance seen by the gate is 



ch 



(11) 
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but, from|2ja), Ctg can also be expressed in terms of capacitances Cg, Csub, and Cg'. 

„ _ Cg {Csub + Cg) 

^tg — ^ ^ I ri ■ ^ ' 

From eqs. ( fTTI ) and ( fT2] i we get the derivative of the channel potential with respect to the gate potential 

(13) 



d(j)ch _ Cg 



Cg + Csub + Cq 

The expression of the sub-threshold slope S then turns out to be 

^=(l + ^^^|^)^ln(10), (14) 

from which it is clear that S is an increasing function of Cq, and therefore a large quantum capacitance severely 
limits device performance. 

|2jb) shows the capacitance-gate voltage characteristics for Vds = 0, 0.1 and 0.25 V obtained by solving the 
Schrodinger equation self-consistently with Poisson equation. In the fully ballistic case the quantum capacitance is 
low for small Vds, indicating a good control of the channel by the gate voltage, but, as soon as Vds increases, hole 
accumulation in the channel occurs and Cq increases, rapidly degrading S (Etb))- In the inelastic case, instead, the 
hole accumulation process is slightly suppressed by inelastic injection from the source but the effect on the quantum 
capacitance is practically negligible. We have observed that for t larger than 10~^ ns, the quantum capacitance 
basically does not change with respect to |2l on the other hand, for r < lO^"* ns Cq decreases with respect to the 
fully ballistic case, but the inelastic process becomes dominant and eq. [14] loses validity. 



IV. Perspectives for device operation 

In order to evaluate the possible performance of the SiC-graphene FET, we have computed the transfer charac- 
teristics by varying three device parameters: drain-source voltage Vds = (ms ^ Mf)/?' donor molar fraction at the 
contacts ao (and therefore parameter A) and oxide thickness tox- We also account for different possible values of 
inelastic time r. First,in [3] we analyze the trend of the transfer characteristics for different t for Vds = 0.25 V, 
tox = 1 nm and an = 6.5 x lO^'* (corresponding to A = 0.01 eV). We observe that for r > 1 ns the transfer 




Fig. 3. Transfer-characteristics for t = 10"^, 10"^, 10"*, 10"^, 10"'' ns, for V^s = 0.25 V an = 6.5 X 10"* and tox = 1 nm. 



characteristics are unaffected and identical to the ballistic case (r oo). Reducing the relaxation time under 1 ns, 
the minimum current increases and the sub-threshold slope remains almost constant since the quantum capacitance 
of the channel does not change. The introduction of inelastic scattering process has mainly two effects in the 
transfer characteristics; one is a gradual change of the current in the sub-threshold region, the other is an increase 
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of saturation current for r < 10^ ns or when inelastic current becomes relevant. In the most favorable case a 
sub-threshold slope of 140 mV/dec can be obtained. 

In|4]we have highlighted the effect of Vds and of the doping level of contacts. As expected the main visible effect 
of increasing Vds is a gradual degradation of the sub-threshold slope, both in the fully ballistic case (|4fa)-(b)) and 
in the case of relaxation time t = 10 ps (|4fc)-(d)), from 84 mV/dec to 202 mV/dec. The reason is simply related to 
the increased accumulation of holes in the channel with increasing Vds, which implies a larger quantum capacitance 
of the channel and therefore a reduced control of the channel potential from the gate voltage. 




Increasing the doping causes an increase of both the maximum current, due to an improved capacity of the 
source to inject electrons, and the minimum current. From |4] we draw the indication that by reducing doping at 
the contacts we improve the current dynamics. As already noted, when the source-drain voltage exceeds the gap of 
the semiconducting channel (Vds > 0.26 V), the characteristics drastically degrade, since band-to-band tunneling 
current becomes comparable with the thermionic current, and hole accumulation in the channel inhibits channel 
control from the gate. 

The increase of oxide thickness tox has mainly two effects, which can be associated to a reduction of the capacitive 
coupling between gate and channel: it increases the sub-threshold slope S (as shown in EJa)), and the opacity 
of tunneling barriers (i.e. a larger A). The former effect is more evident for Vds ~ 0.25 V, where the quantum 
capacitance is larger, instead S is almost constant at about 75 mV/dec for Vds = 0.1 V and r < lO^"' ns, instead 
for smaller t, S ^ 100 mV/dec. |5jb) represents the /on/^off ratio as a function of for Vds = 0.1 V and 
Vds = 0.25 V, calculated for a gate voltage range AVg — 0.25 V for different values of r. Larger values of the 
lon/IoS ratio are observed for r = 1 ns. 

From our analysis of transfer characteristics, evaluated by varying three device parameters as Vds, ctD and tox, we 
stress the important result that for small Vds and doping level (Ufa)) current is modulated by more than 4 orders of 
magnitude. 

We have to stress also the main limitation of graphene on SiC: the energy gap of 0.26 eV coupled to a low 
effective mass results in a high band-to-band tunneling current Vds > 0.26 V and so in an increase of the minimum 
current achievable. This limitation on Vds affects the perspectives for digital circuit operation: in that case we need 
Vds — AVg and equal to the supply voltage. Even for optimized device parameters (tox — 1 nm, = 4 x 10^"^), 
and a supply voltage of 0.25 V, we obtain an /on/^off ratio of 50, as can been seen from |3b). 
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Fig. 5. a) Sub-threshold slope for V^s = 0.1 V and Vda = 0.25 V calculated with different relaxation time r. b) /on/^off ratio for Vjs = 0.1 V 
and Vds = 0.25 V calculated with different relaxation time r. 



V. Conclusion 

In this work we have investigated the performance of field-effect transistors based on epitaxial graphene on a SiC 
substrate with an analytical model. We have shown that, for small Vds and doping level, current is modulated by 
more than four orders of magnitude: this is a main improvement with respect to other graphene-based devices [20], 
[21], [16], [15]. Comparable results can be obtained only with carbon nanotubes or graphene nanoribbons, but only 
with post-selection of devices after fabrication (for proper chirality and/or width). In the case of graphene on SiC, 
lithography and device patterning are certainly not prohibitive. A steep subthreshold behavior (5* = 67 mV/decade) 
can be obtained for small Vds = 0.1 V, when the accumulation of holes in the channel is inhibited, and a larger 
current ratio, in excess of 10"^, can be obtained for a gate voltage window of 0.25 V. For digital applications, the 
limiting factor is represented by the small voltage drop applicable to the channel, being limited by the energy gap 
(0.26 eV) of the semiconducting material. With optimized device parameters we have obtained a sub-threshold slope 
of « 140 mV/decade and an lon/IoS equal to 50, with a supply voltage of 0.25 V and r = 1 ns. This falls short 
of requirements of the International Technology Roadmap for Semiconductors, which requires lon/IoS ~ 10^ [22]. 
Finally, we believe that graphene on SiC is very promising as a channel material for FETs, and much attention has 
to be put on mechanisms capable to suppress hole injection also at larger Vds, that would allow to improve the 
subthreshold swing and obtain a good /on/^off also with a small applied voltage, and on its use in tunnel FETs, 
where its low gap and low effective mass can be turned into an advantage. 
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